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! ABSTRACT 

We present results from a study of the X-ray cluster population that forms within 
the CLEF cosmological hydrodynamics simulation, a large A^-body/SPH simulation 
of the ACDM cosmology with radiative cooling, star formation and feedback. With 
nearly one hundred {kT > 2 keV) clusters at z = and sixty at z = 1, our sample is 
] one of the largest ever drawn from a single simulation and allows us to study variations 

within the X-ray cluster population both at low and high redshift. The scaled projected 
temperature and entropy profiles at z = are in good agreement with recent high- 
quality observations of cool core clusters, suggesting that the simulation grossly follows 
Q ■ the processes that structure the intracluster medium (ICM) in these objects. Cool cores 

^ I are a ubiquitous phenomenon in the simulation at low and high redshift, regardless of a 

C/3 . cluster's dynamical state. This is at odds with the observations and so suggests there is 

' still a heating mechanism missing from the simulation. The fraction of irregular (major 

merger) systems, based on an observable measure of substructure within X-ray surface- 
brightness maps, increases with redshift, but always constitutes a minority population 
' within the simulation. Using a simple, observable measure of the concentration of the 



H I ICM, which correlates with the apparent mass deposition rate in the cluster core, 

we find a large dispersion within regular clusters at low redshift, but this diminishes 
at higher redshift, where strong cooling-flow systems are absent in our simulation. 
Consequently, our results predict that the normalisation and scatter of the luminosity- 
temperature relation should decrease with redshift; if such behaviour turns out to be a 
correct representation of X-ray cluster evolution, it will have significant consequences 
for the number of clusters found at high redshift in X-ray flux-limited surveys. 
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1 INTRODUCTION 

Clusters of galaxies, occurring at low redshifts (z < 2), 
are interesting cosmological objects as they offer a 
powerful yet independent approach from other meth- 
ods (such as the cosmic microwave background, at 
z ~ 1000) for constraining cosmological parameters. 
While several cosmological applications of clusters exist, 
a particularly appealing method, because of its sim- 
plicity, is to measure the va riation in the cluster mass 
function with redshift (e.g. iBlanchard &: BartlettI 1 19981 : 
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lEke et al.l 1 19981 ). Since mass is known to be tightly 
correlated with X-ray observables, particularly tem- 
perature (e.g. Finoguenov. Reiprich fc Bohringer 2001 



Arnaud. Pointecouteau fc PrattI l2005l : IVikhlinin et al.l 
2006al ). it is straightforward in principle to convert between 



the two quantities. 

Theoretically, such cluster scaling relations were pre- 
dicted to exi st, esse n tially as a manifestation of the virial 
theorem, by iKaiseij (|l986l ). In the so-called gravitational- 
heating scenario, the intracluster medium (ICM) was heated 
by the gravitational collapse and subsequent virialisation of 
the cluster. X-ray observations of clu sters confirmed the ex- 
istence of these scaling relations (e.g. lEdge fc Stewartlll99ll : 
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iFabian et alll 19941 for the X-ray luminosity-temperature re- 
lation) though they revealed two complications. Firstly, the 
slope of the observed X-ray luminosity-temperature rela- 
tion (and to a lesser extent, the mass-temperature relation) 
is steeper than pre dicted from gravitational heat i ng alo ne, 
an effect shown by iPonman. Cannon fc Navarrol (|l999l ') to 
be due to an excess of entropy in the cores of clusters, and 
more so in groups. A lot of theoretical effort has gone into 
understanding the origin of the excess entropy (see lVoitlbOOSi 
for a recent review). Secondly, there is an intrinsic scatter in 
the scaling relations, which is particularly large for the low- 
redshift luminosity-temperat ure relation, due to the large 
variations in core luminosity (jPabian et al.|[l993 ). For cos- 
mological studies with clusters, an accurate statistical de- 
scription of the cluster population is warranted, as only then 
can robust cluster survey selection functions be constructed. 
From a theoretical stand-point, cluster scaling relations of- 
fer an additional, exciting prospect; the amount with which 
these relations evolve with redshift ought to reveal informa- 
tion on the nature of no n-gravitational processes and clust er 
astrophysics in general (|Muanwong. Kav fc ThomaslbOO^ ). 

The intrinsic scatter in cluster scaling relations can 
at least partly be attributed to gravitational processes, 
as cl usters themselves liv e in different environments 
(e.g. ISchuecker et all l200ll ). although non-gravitational 
processes, such as radiative cooling and heating from 
galaxies, must also p l ay a role (e. g. Pearce et al.l I2OO0I : 
iMcCarthv et al.l l2004l : iBalogh et aP I2OO6I ). Cosiiiological 
simulations of the cluster population are the most ac- 
curate method by which to characterise the statisti- 
cal properties of clusters, as they include an accurate 
treatment of the non-linear gravitational dynamics and 
merging processes, as well as allowing non-gravitational 
physics to be incorporated self-consistently. Early at- 
tempts focused on the simplest model for the gas, a non- 
radiative ICM, which was successfully shown to repro- 
duce the simple, self-similar, sca lings expected from the 
gravi t ational-heating model (e.g. iNavarro. Frenk fc White 



1995; 'Evrard, Metzler & NavarrJ 'l996'; 'Brvan & Norman' 
1998; Eke, Navarro & Frcnk 1998; Mu anwong et al. ,200 2). 

Additional non-gravitational processes have also been 
studied within simulations, and various mechanisms 
have been proposed to e xplain the similarity breaking , 
such as preheating (e . g. [ Navarro, Frenk fc Whit3 1 19951 : 



iBialek. Evrard fc Mohr 



2001: 



B organi et al 

five c ooling fe.g. [ Pearce et al.l I2OO O: Mu anwong et al 



20021). radia 



2001 . l2002l: iDave. Katz fc WeinbergI 120021 : iMotl et al 



2004: iKravtsov Nagai fc VikhhninI 120051 ') or both (e.g 



Muanwong et alj 120021 )7" Recently, attention has shifted 



to more realistic models which attempt to directly cou- 
ple feedback (local hea ting from galaxi e s) with cooling 
and s tar formation (e. g. Valdarnini zOCXj; 'To rnatore et al.l 



20031: iKav. Thomas fc Thcuns 2003; Borga ni et al.l |2004 



Kavll2004l : iKav et al.ll2004l : lEttori et al.ll2004l ) 



Together with progress in the development of these non- 
gravitational models, the advance in both simulation codes 
and computer hardware is now allowing larger simulations 
with reasonable resolution to be performed. We are now be- 
ginning to resolve sufficient numbers of clusters to start mak- 
ing quantitative predictions at all appropriate redshifts for 
the cluster population. The CLEF-SSH (CLuster Evolution 
and Formation in Supercomputer Simulations with Hydro- 



dynamics) collaboration has been set up to take advantage 
of this new era in numerical modelling, by performing large 
simulations of the cluster population. Our first simulation, 
known as the CLEF simulation, is a large (A'^ = 2 x 428^) 
A''-body/SPH simulation of the ACDM cosmology, within a 
200/i~^Mpc box, and includes a model for radiative cool- 
ing and energy feedback from galaxie s. This simulat i on is a 
similar size to the one performed bv lBorgani et all 
but uses a different feedback model. In iKav et al.l 



2004), 



20051 ). 



hereafter Paper I, we presented a small selection of results 
a,t z — from the CLEF simulation. For this paper, we have 
performed a more detailed analysis of the same cluster popu- 
lation, and present results for a range of redshifts from z = 
to 1, focusing on the effects of dynamical activity and the 
strength of cooling cores. A companion paper (da Silva et al., 
in preparation) presents results for the Sunyaev-Zel'dovich 
effect properties of the CLEF cluster population. 

The rest of this paper is organised as follows. In Sec- 
tion [5] we summarise details of the CLEF simulation and 
detail our method for creating cluster catalogues, maps and 
profiles. The internal structure of clusters, and how it de- 
pends on dynamical regularity and the properties of the 
core, is the focus of Section|3] Section [4] then draws on these 
results to investigate the evolution of key cluster scaling re- 
lations with redshift. We discuss our results in Section [S] and 
summarise our conclusions in Section [6] 



2 THE CLEF SIMULATION 

The CLEF simulation is a large (A^ = 2 x 428^ particles 
within a 200 ^~^Mpc comoving box) cosmological simulation 
of structure formation, incorporating both dark matter and 
gas. Below we describe the procedure used to generate the 
simulation data and how the clusters were identified within 
these data to create X-ray temperature- limited samples from 
redshifts z — to z — 1. 



2.1 Simulation details 

For the cosmological model, we adopted the spatially-flat 
ACDM cosmology, setting the following values for cosmo- 
logical parameters: matter density parameter, f2m ~ 0.3; 
cosmological constant, Qa ~ A/3Hq — 0.7; baryon den- 
sity parameter, Qh = 0.0486; Hubble constant, h = 
-ffo/lOOkms^^Mpc^^ — 0.7; primordial power spectrum in- 
dex, n = 1 and power spectrum normalisation, as — 0.9. 
These values were chosen to be consistent with the results 
from WMAP first year data (|Spergel et al.ll2003l ). 

Initial conditions were generated for a cube of comoving 
length 200/!-~^Mpc at redshift, z = 49. The cube was popu- 
lated with two interleaving grids of 428'^ particles, one grid 
representing the dark matter and one representing the gas; 
the particle masses were thus set to rridark = 7.1 x 10® h~^MQ 
and rrigas = 1.4 x 10® h~^MQ for the dark matter and gas re- 
spectively. Initial particle displacements and velocities were 
then computed fro m a transfer function genera ted using the 
CMBFAST code jSeliak fc Zaldarriagal Tl996l ) . The initial 
temperature of the gas was set to T = 100 K, significantly 
lower than the range of temperatures typical of overdense 
structures resolved by the simulation. 
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The initial conditions were evolved to z = usinj 



n g a ver - 

sion of the GADGET2 iV-body/SPH code (jSpringefeool ). 
modified to include additional physical processes (radiative 
cooling, star formation and energy feedback; see below). 
Gravitational forces were calculated using the Particle-Mesh 
algorithm on large scales (using a 512^ FFT) and the hierar- 
chical tree method on small scales. The (equivalent) Plum- 
mer softening length was set to e = 20/i~^kpc, fixed in co- 
moving co-ordinates, thus softening the Newtonian force law 
below a comoving separation, a;min ~ 2.8 e = 56ft~^kpc. 

Gas particles were additionally subjected to adia- 
batic forces, and an artificial viscosity where the fiow was 
convergent, using the entro py-conserving version of SPH 
llSpringel fc Hernguistl 120021 '). the default method in GAD- 
GET2. Additionally, we allowed gas particles with T > lO^K 
to cool radiativel y, using th e isochoric cooling approxima- 
tion suggested by iThomas fc Couchman ( 199 ^). Tabulate d 
cooling rates were taken from Sut herland fc Dopital (|l993h , 
assuming an optically thin Z — 0.3 Zq plasma in coUisional 
ionisation equilibrium (a good appro ximation to the i ntra- 
cluster medium out to at least z = l: lTozzi et al.|[20o3 ). 



2.2 Feedback 

We have also attempted to follow crudely the large-scale 
effects of galactic outfiows (feedback) in the simulation, to 
regulate the cooling rate and to inject non-gravitational en- 
er gy into the gas. We adopted the Strong Feedback model 
of lKavl l|2004l ). hereafter K2004a, as it was shown there 
to approximately reproduce the observed excess entropy in 
groups/clusters, both at small and large radii (see also Pa- 
per I). We only give a brief summary of the model details 
here. 

First of all, cooled gas is identified with overdensity 
5 > 100, hydrogen density nn > 10"'^ cm~'^, and temper- 
ature T < 12, OOOK. For each cooled gas particle, a random 
number, r, is drawn from the unit interval and the gas is 
reheated if r < /heat, where /heat = 0.1 is the reheated mass 
fraction parameter. Reheated gas is given a fixed amount 
of entropyQ Sheat ~ lOOOkeVcm^, corresponding to a mini- 
mum thermal energy of ~ 17keV at the star- formation den- 
sity threshold. Such a high thermal energy (compared with 
typical cluster virial temperatures) means that the reheated 
gas is supersonic and is thus distributed through viscous in- 
teractions and shocks in the ICM. T his not only regula tes 
the star-formation rate in the cluster l|Balogh et al.|[200ll ). it 
also prevents sig nificant build-up of low-entropy material i n 
the cluster core (|Kav. Thomas fc Theunsll2003l : lKavll2004 ). 
Our model could thus be perceived as a crude representa- 
tion of local accretion-triggered heating by stars and active 
galactic nuclei (although feedback from the latter does not 
necessarily have to follow the star-formation rate, as is done 
here) . 

Gas particles that are not reheated are instead con- 
verted to coUisionless star particles. Although the model 
does not treat star formation (which occurs in regions with 



^ We define entropy as 5 = kT{p/^mii)^ where 7 = 5/3 
the ratio of specific heats for a monatomic ideal gas and /imn 
0.6 is the mean atomic weight of a fully-ionised plasma. 



much higher gas densities, nn ^ 0.1) accurately, this pre- 
mature removal of low-pressure material from the gas phase 
saves computational effort as these particles generally have 
the shortest timesteps. Furthermore, it helps to alleviate the 
difficulty that standard SPH has in resolving the sharp in- 
terface between hot and cold phases l|Pearce et al.ll2()()0l ). 



2.3 Cluster identification 

The CLEF simulation produced a total of 72 snapshots of 
the particle data, at time intervals optimised for producing 
mock lightcones (da Silva et al., in preparation). Only the 
25 lowest-redshift snapshots are used in this paper, rang- 
ing from 2; = to 1; at higher redshift the number of clus- 
ters becomes prohibitively small. We used these snapshots to 
produce cluster catalogues (mass, radius and various other 
properties), maps and profiles. 

Catalogues we re generated using a si milar procedure 
to that adopted bvlM uanwong et al.l (|2002l '). Briefly, groups 
of dark matter particles we re identified using the friends- 
of- friends (FoF) algorithm (jPavis et al.l 119851 ) . setting the 
dimensionless linking- length to b = 0.1. Spheres were then 
grown around the particle in each group with the most nega- 
tive gravitational potential, until the enclosed mass equalled 
a critical value 



(1) 



where A is the density contrast, pcr(z) = {SHq /8'kG)E{z)^ 
is the critical d ensity and E(z)^ = Om( l + z)^ -I- 1 — fim for 
a flat universe ijBrvan fc Normanlll998h . Cluster catalogues 
for a variety of density contrasts were constructed, including 
th e virial value (A 100 at z — 0), taken from equation (6) 
in iBrvan fc Normaiil (|l998h . The virial radius was used to 
find overlapping pairs, and the least massive cluster in each 
pair was discarded from the catalogues. For nearly all of 
the results presented in this paper, we use a catalogue with 
A — 500, as this is the smallest density contrast typically 
accessible to current X-ray observations. 

We initially selected all clusters with at least 3000 par- 
ticles within i?5oo, corresponding to a lower mass limit of 
Afsoo = 2.5 X 10^^ /i"^M0. This limit is low enough that our 
temperature-selected sample (below) is comfortably a com- 
plete subset of this sample at all redshifts studied. At z = 
we have 641 clusters in our mass-limited sample, decreasing 
to 191 clusters at z = 1. This is comp arable to the number s 
found in the simulation performed bv lBorgani et al" I (|2004h . 
who also used GADGET2 but with a different prescription 
for cooling, star formation and feedback than used here. 



2.4 Spectroscopic-lilce temperature 

Observational samples of X-ray clusters are usually limited 
in flux or temperature, where the latter is measured by fit- 
ting an isothermal plasma model to the observed spectrum 
of the cluster X-ray emission. Theoretical models of X-ray 
clusters commonly use an emission- weighted temperature to 
estimate the spectral temperature of a cluster. For particle- 
based simulations, this is done using the formula 



T : 



(2) 
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Figure 1. X-ray emission-weighted temperature, Tew, plotted 
against spectroscopic-like temperature, Tgi, for clusters with 
fcTgi > 2 keV. The solid line corresponds to Tew = T^i and the 
dashed line is the best-fit straight line to our data. The dotted 
line is the best fit to the clusters studied by Rasia et al. (2005). 



2.5 Cluster maps and projected profiles 

Cluster maps w ere produced usinR a similar p rocedure to 
that discussed in lOnuor;i7 Kav fc Thoma3(|2003l ): m essence, 
values at each pixel are the sum of smoothed contributions 
from particles, using the same spline kernel as used by the 
GADGET2 code. Centred on each cluster, only particles 
within a cube of half-length, I — 4i?5oo, were considered, 
i.e. out to approximately twice the virial radius in each or- 
thogonal direction. For this paper, we computed bolometric 
surface brightness (although the emission is predominantly 
thermal bremsstrahlung in the X-ray) and spectroscopic-like 
temperature maps. Projected temperature and azimuthally- 
averaged surface-brightness profiles were also computed, av- 
eraging particles within cylindrical shells, centred on the 
pixel with the highest surface brightness. 

Fig. [2] illustrates spectroscopic-like temperature maps 
of the 5 most massive clusters eac h at ^ = and z = 1, out 
to a r adius, -R500 . As w a s foun d bv lOnuora. Kav fc Thomas! 
('2003) and iMotl et al.l (|2004l ). there is a large amount of 
temperature structure within each cluster, particularly cold 
spots due to cool, low-entropy gas trapped within infalling 
sub-clusters. The intensity scale is defined by the minimum 
and maximum temperature; the dynamic range is typically 
an order of magnitude (Tmax/Tmin), with maximum temper- 
atures being around twice that of the mean. 



where Wi = miniA{Ti, Z) is the weight given to each hot 
(Ti > 10^ K) gas particle i, rrn its mass, n,; its den- 
sity, Ti its temperature and K{Ti,Z) the cooling func- 
tion, usually for emis s ion within the X-ray ene rgy band 
(e.g. iMuanwong et al.l l2002l : iBorgani et al] 120041 ) . In the 
bremsstrahlung regime A oc T^^'^ and so the hottest, densest 
part icles are given t h e mos t weight. 

iMazzotta et al.l l|2004l ) applied the same method used 
by observers to measure the spectroscopic temperature of 
simulated clusters, and found that it was always lower than 
Tew. For the bremsstrahlung regime {kT > 2 keV), they sug- 
gested a more accurate measure, known as the spectroscopic- 
like temperature, with weight Wi = miniT~'' . This esti- 
mator gives the coldest, densest particles more weight. We 
adopt this estimator in this paper (summing over gas par- 
ticles within -R500 with kTi > 0.5 keV, the typical lower en- 
ergy limit of an X-ray band) to create a temperature- limited 
(fcTsi > 2 keV) sample of clusters at all redshifts. This re- 
duces the number of clusters to 95 at z = 0, decreasing to 
57 at z = 1. While this is one of the largest temperature- 
selected cluster sample drawn from a single simulation, we 
note that the dynamic range is still quite small. Nearly all 
clusters have Tsi ~ 2 — 4 keV at all redshifts (the median 
temperature stays approximately constant with redshift at 
Tbi ~ 2.5 keV) and our hottest cluster at 2; = has T^i = 7.3 
keV. 

Fig. [1] compares Tgi to Tew (for a 0.5 — 10 keV band) 
for our temperature- limited sample at 2: = 0. As was found 
by Rasia ct al. ( 20051) , whose sample mainly consisted of the 
[Borgani ct al. ( 2004 ) clusters, Tsi and Tew differ by as much 
as 20 per cent. Rasia et al. found fcTsi = 0.7fcrew + 0.3, 
whereas we find kT^i = O.SfcTew + 0.1, sim ilar to, although 
slight ly steeper than, their result (see also iKawahara et al.l 
I2OO7I ). 



3 CLUSTER STRUCTURE 

In this section we present the structural properties of the 
CLEF clusters, comparing to observational data where ap- 
propriate. 



3.1 X-ray temperature bias 

As discussed in the previous section, the X-ray temperature 
of the ICM is biased to regions of high density. Cooling and 
heating processes are generally most efficient there, so the 
X-ray temperature of a cluster is not necessarily an accurate 
measure of the depth of the underlying gravitational poten- 
tial well, even if the system is virialised and approximately 
in hydrostatic equilibrium. 

We investigate any such temperature bias in our simula- 
tion by comparing T^i to the dynamical temperature, Tdyn, 
through the standard quantity, /Jspcc ~ Tdyn/Tai. The dy- 
namical temperature. 



fcTc 



dyn 



V. rrnkTi+aY,. h 



(3) 



where ct — (2/3)/imH ~ 6.7 x 10~^^ g, assuming the ratio of 
specific heats for a monatomic ideal gas, 7 = 5/3, and the 
mean atomic weight of a zero metallicity gas, /irriH = 10"^* 
g. The first sum in the numerator runs over gas particles and 
the second sum over all particles, of mass rrii, temperature 
Ti and speed Vi in the centre of momentum frame of the 
cluster. 

Fig. [3] illustrates /Jgpcc for each cluster at 2 = 1 and 
2 = 0, versus its scaled mass, _E(2) A/500 • Only clusters with 
E[z)M5oo > lO"/i-'M0 are selected, producing similar 
numbers to our temperature-selected samples at both red- 
shifts. There is a clear positive correlation between /Jspcc and 
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Figure 2. Spectroscopic-like temperature maps of the 5 most massive clusters (in order of decreasing mass, left to right) at z = 1 (top 
panels) and 2 = (bottom panels). The spectroscopic-like temperature is given in each panel. Surface- brightness contours (normalised to 
the maximum value) are overlaid; adjacent contours correspond to a difference in surface brightness of a factor of 4. Images are centred 
on the maximum surface-brightness pixel and are shown out to -Rsoo- 
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Figure 3. Ratio of dynamical temperature to spectroscopic-like 
temperature, /Jspoc = Tdyn/Tgi, versus scaled mass, E(2)M5oo, at 
z = 1 and 2 = 0. The solid line is the best-fit relation to the data, 
while the dot-dashed (dashed) lines are best-fit relations when Tgi 
is replaced by the hot gas mass-weighted temperature excluding 
(including) bulk kinetic motions. 



i?(2) A/500, with /3spec > 1 for most clusters (i.e. Tgi < Tdyn) 
at low and high redshift. The spectroscopic-like temperature 
is a biased tracer of the gravitational potential for three rea- 
sons. Firstly, cool dense gas is weighted more than less dense 
material, as discussed in the previous section. This effect can 
be seen in the figure by comparing jSspcc to the best-fit rela- 
tion when Tai is replaced by the hot gas mass-weighted tem- 
perature (dot-dashed line). Secondly, some of the energy of 
the gas is in macroscopic kinetic energy, as can be deduced 
from comparing the dot-dashed to the dashed line, where 
in the latter case, T^i is replaced by the temperature when 
equation (|3} is applied to only the hot gas. Finally, feedback 
heats the gas, particularly in low mass clusters (the dashed 
line shows that /3spcc < 1 for most clusters, i.e. the gas has 
more specific energy than the dark matter). 



3.2 Baryon fractions 

We also examine the segregation of baryonic mass into 
gaseous (ICM) and galactic (coUisionless) components 
within each cluster. Gas and baryon fraction profiles for this 
model have already been studied bv lKav et al.l (|2004l ). here- 
after K2004b, who showed that the baryon fraction profiles 
were in good agreement with observations but that too much 
of the gas had turned into stars (the normalisation of the gas 
fraction profile is as low as 50 per cent of the observed pro- 
file). Here we examine the behaviour of the baryon/gas/star 
fractions with temperature and redshift. 

Fig. m illustrates baryon fractions normalised to the 
global value, Yh = /b/(fib/nm) within R500, for each cluster 
at z = 1 and z — 0. Most of our clusters have feTsi < 5 keV, 
where there is a strong trend in increasing baryon fraction 
with temperature, as feedback can heat and expel more gas 
in smaller clusters. At high temperature, a few systems at 
z = are consistent with the mean value (~ 0.9) found by 
K2004b for their non-radiative clusters. Overall, the mean 
baryon fraction increases by 8 per cent between the two red- 
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Figure 4. Baryon (squares), gas (triangles) and star (diamonds) 
fractions versus spectroscopic-like temperature at z = 1 and 
z = 0. Horizontal dashed lines illustr ate mean values. The solid 
line is the average value measured bv lKav et al. 1 ||2004|') for their 
non-radiative clusters. 



Figure 5. Substructure statistic, Sx, versus spectroscopic-like 
temperature for clusters at z = I and 2 = 0. Triangles illustrate 
regular clusters with Sx ^0.1 and squares irregular clusters with 
Sx > 0.1. The solid horizontal line is the median Sx and the 
dashed lines the 10 and 90 percentiles. 



shifts, from 0.71 at 2; = 1 to 0.79 at 2; = 0. Similarly, the 
mean hot gas fracti on increases f r om .42 to 0.49 over the 
same redshift range. lEttori et all l|2004h also found the gas 
fractions to weakly decrease with redshift in their simulated 
clusters, albeit with higher values than found here. 

The star fraction is a very weak function of both tem- 
perature and redshift, with a mean value of 0.28 at 2 = and 
0.29 at z — 1. Just under 40 per cent of the baryons within 
R500 have condensed and formed stars in our simulation, at 
all redshifts; a value that only decreases to about 30 per cent 
at the virial radius. Observations indicate a value of about 
10 to 15 per cent, significan tly lower than in our clusters 
l|Lin. Mohr fc Stanford! [20031 ). Thus, as found by K2004b, 
our feedback model has not been effective enough at limit- 
ing th«_overcooling_of baryons in clusters, as was also found 
bv lEttori et al.l l|2006l) . However, the global star fraction is 
only 13 per cent at z = (and 9 per cent at 2 = 1), just 
sligh tly larger than the observed value of 5 to 10 per cent 
(e.g. iBalogh et allbOOll ). 

3.3 Regularity 

Hierarchical models of structure formation predict that sub- 
structures in clusters should be commonplace, as clusters 
are the latest result of a series of mergers of smaller sys- 
tems, and have dynamical times (~ 1 Gyr) that are a sig- 
nificant fraction of the age of the Universe. Indeed, sub- 



structure is frequently observed in clusters and dynami- 
cal activity ha s been quantified using various techniques 
(e.g.lJones fc F orman 1992; Mohr, Fabricant & GcUcr 199^; 
Buote fc T sai 1995, 1996; Cr one. Evrard fc Richstonelll996l : 
Schuecker et al...2001 : Jchcma et al.ll2005h . 



We use a simple measure of substructure in our clus- 
ter surface-brightness maps, using the centroid-shift method 
simila r to that suggested by .Mohr. Fabricant fc Gelleil 
(|l993h 



Sx = 



|Rz 



Us 



(4) 



where Re, max is the position of the pixel with maximum 
surface brightness (taken to be the centre of the cluster) 
and R s.cen is the surface-brightness centroid. iThomas et al.l 
(ll99St ) used a similar method, based on the 3D total mass 
distribution, and showed that this did as well, or better, 
than other more sophisticated measures of substructure in 
simulated clusters. 

Fig. [5] illustrates Sx versus T^i for our temperature- 
limited sample of clusters ai z — 1 (top panel) and 2 = 
(bottom panel). It is evident that the range of Sx values at 
fixed temperature is large: the 10 and 90 percentiles of each 
distribution, shown as horizontal dashed lines, vary from 
~ 0.01 to ~ 0.15. Inspection of surface-brightness maps (see 
Fig. [6]) reveals that clusters with the largest Sx appear dy- 
namically disturbed and are therefore undergoing a major 



X-ray clusters in the CLEF simulation 7 



kT=2,6keV S<=0,00 


kT=2,2keV S<=0,00 




kT=2,4keV S<=0,01 


kT=3,4keV S,=0,01 






kT=Z4keV^__Sj=0,35 




(^} 





Figure 6. Bolometric surface-brightness maps of the 5 clusters with the lowest (top panels) and 5 with the highest (bottom panels) 
substructure statistic (5x) at z = 0. In this case, clusters are centred on the surface-brightness centroid. 



merger. We choose to divide our sample into irregular clus- 
ters with Sx > 0.1 and regular clusters otherwise. We note 
that this division is somewhat arbitrary and only serves to 
provide us with a means to compare the most disturbed clus- 
ters at each redshift to the rest of the sample. The longer tail 
in the Sx distribution to high values at z — exacerbates 
the difference between the two sub-populations relative to 
those at z = 1 (the length of the tail itself changes from 
redshift to redshift). 

There is no significant trend in Sx with temperature, 
within the limited dynamic range of our sample. However, 
there is a trend in Sx with redshift: the median value at 
z = 1 is almost a factor of 2 higher than at 2 = 0. In other 
words, clusters tend to be less regular at higher redshift. 

The increase in dynamical activity with redshift in our 
simulated cluster population is qualitat ively consistent with 
the recent result of lJeltema et al.l (|2005|), w ho used the more 
complex power ratios l|Buote fc Tsail 199^ to measure dy- 
namical activity in a sample of low- and high-redshift clus- 
ters observed with Chandra. 



3.4 Temperature and surface-brightness profiles 

Surface-brightness and projected temperature profiles are 
now regularly observed for low-redshift clusters with XMM- 
Newt on and Chandra (e .g. Arnaud. Pointecouteau fc Pratt! 



20061 ; 



20051: Ivikhhnin et al.l |2005|: 



Zhang et al.l 

observable quantities, as 
information can be extracted 



Piffaretti et all l2005l : 
Pratt et all 120071 ). These are key 
3D density and temperature 
from these measurements 



through deconvolution techniques. This allows the thermo- 
dynamics of the ICM to be studied, as well as the total 
mass distribution to be calculated (assuming hydrostatic 
equilibrium) . 

In Fig. [7] we present scaled projected spectroscopic- 
like temperature profiles for our regular and irregular 
clusters at z = 1 and 2 = 0. As stated previously, each 
cluster (including irregular objects, where all emission 
is included in our analysis) is centred on the pixel with 



maximum surface brightness. For consistency with the 
observational data, each temperature profile is nor- 
malised to the average spectroscopic-like temperature, 
(Tsi), between projected radii of 50/i~^kpc and i?5oo- 
Projected radii are then r e-scaled to R iso using the 
formula Riso = l.gSi/fc {T^i} /lOkeV /E{z) ft~^Mpc, 
originally derived from numerical simulations by 
lEvrard. Metzler fc Navarrol l|l996l ). 

As was found by K2004b, the median profile rises 
sharply from the centre outwards, peaks at ~ O.l-Rigo, then 
gradually declines at larger radii. The inner rise, where the 
density is largest, is due to radiative cooling of the gas, while 
the outer decline is a generic prediction of the ACDM model 
(e.g. Eke, Navarro fc Fronk 1998). It is interesting to note 
that the shape of the profile for irregular clusters is flatter 
than for the regular majority, beyond the peak. This is due 
to the second, infalling, object, which compresses and heats 
the gas. We also note that the temperature profiles ai z — 1 
are very similar to those at 2 = 0, and so a cool core is 
established in the cluster early on. 

IVikhlinin et al.l (|2005l ) recently determined the pro- 
jected temperature profile for a sample of 11 low-redshift 
cool core clusters observed with Chandra. The shape of their 
profile is very similar, albeit slightly steeper at large radii, 
to that of our regular clusters; a rough fit, as sup plied by the 
authors, is shown in Fig. [7] as thick solid lines. iPratt et all 
(2007) performed a similar study with XMM-Newton, for 
a sample of 15 clusters (including non-cool core systems); 
their result is shown in the fi gure as the shaded region. In- 
terestingly, iPratt et all (|2007l ) find a similar decline at large 
radius to our regular clusters but the temperature does not 
drop as sharply in the centre (even for those clusters with 
coolest cores). 

The presence of cool cores at both low and high redshift 
in our simulation is i n qua litative agreement with the find- 
ings of iBauer et al.l (|2005l ). who measured central cooling 
times for a sample of 2 = 0.15 — 0.4 clusters observed with 
Chandra and found their distribution to be very similar to 
that for a local sample. 
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Figure 7. Scaled projected spectroscopic-like temperature pro- 
files at z = 1 and 2 = 0. Solid curves are median and 10/90 
percentiles for regular clusters, and dashed curves for irregular 
clusters. The vertical dashed line illustrates the median scaled 
softening radius (i.e. where the gravitational force becomes softer 
than Newtonian). The thick solid lines (with zero and negative 
gradient) are fits to the average o bserved temperature p rofile of 
cool core clusters, as measured bv lVikhlinin et al.l ||2005|) : the in- 
ner line is a rough fit to their data to illustrate the cool core. 
The shaded region encloses the mean and la standard deviation 
te mperature p r ofile f or a representative sample of nearby clusters 
bv lPratt et"all ||2007|) . 



Bolometric surface-brightness profiles are presented in 
Fig. [8l At both redshifts, it is clear that there is a larger 
dispersion between clusters in the core than at the outskirts, 
particularly at « = 0. The irregular clusters have fiatter 
profiles than the regular clusters and a bump can be seen at 
large radius, due to the core of the second object. 

We also calculate density and temperature gradients 
for our clusters, as is needed for cl uster mass estimates 
(Section [33. Following Ivikhlinin et~a l. (2006a), we define 
PcS = -(l/3)d Inp/dlnr and A = -(1/3) dluT/dlnr 
to represent 3D density and temperature gradients respec- 
tively. Fig. |9] shows these values for our clusters at i^soo , 
plotted against temperature. Resu lts at z = are o verplot - 
ted with the Chandra data from IVikhhnin et~aLl (l2006al ). 
In general, the agreement between our results and the ob- 
servations is very good; median values are /3eff = 0.76 and 
Pt = 0.19 respectively. At z = 1 the median values change 
very little (0.73 and 0.18). 
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Figure 8. Bolometric surface-brightness profiles at 2 = 1 and 
2 = 0. Again, solid curves are median and 10/90 percentile val- 
ues for regular clusters, and dashed curves for irregular clus- 
ters. The vertical dashed line marks the median force resolution, 
(2.8e/i?500>. 
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Figure 9. Effective slopes of gas density (/3eff ) and temperature 
(/3t) profiles at i?500i versus core-excised spectroscopic-like tem- 
perature, for clusters at 2; = 1 and z = 0. Triangles are regular 
clusters and squares irregular clu sters. Data points with error bars 
are from I Vikhlinin et al'l (Hooli). 
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Figure 10. Spectroscopic-like temperature maps of the 5 clusters with the lowest (top panels) and highest (bottom panels) /x values at 
2 = 0. Circles mark the two inner radii (0.1 and 0.3 -R500) where /x is measured. X-ray concentrations (/l) are also given, which tend 
to be anti-correlatcd with /t for regular clusters. 



3.5 Core structure parameters 

Two generic features of our simulated clusters at 2 = is 
that the majority have cool cores (as shown in Fig. [7]) and 
exhibit a large dispersion in core surface brightness (Fig.lS)). 
At 2 = 1, the clusters also tend to have cool cores but a 
smaller range in core surface brightness is seen. To quantify 
this behaviour further, we define two simple core structure 
parameters that are readily observable. 

We first define three projected radii, \R\,Ri,Ri\ = 
[0.1, 0.3, 1.0]7?5oo. The first approximately defines the radius 
where the profile stops rising; the first and second approxi- 
mately define the (maximum) temperature plateau, and the 
third is the outer radius of the cluster. The first parameter 
is then 



/t 



T,i(< Ri) 



r,i(i?i ^Ri^R2)' 



(5) 



which measures the ratio of the core to the maximum pro- 
jected spectroscopic-like temperature of the cluster. Clusters 
with the coolest cores have the lowest /t values. This can 
be seen clearly in Fig. 1101 where temperature maps of the 5 
clusters with the lowest and 5 with the highest /t values are 
shown. All but two clusters in our sample, including irregu- 
lar systems, have /t < 1 because of their cool cores; at 2 = 1 
the situation is similar, where only 6 clusters (~ 10 per cent 
of the sample) have /t > 1 (the median /t increases gradu- 
ally with redshift). We will return to this point in Section [S] 
The second structure parameter is 



Lhoi{< Ri) 

ibol(< R3)' 



(6) 



which measures the fraction of bolometric luminosity ema- 
nating from the core; we label this the X-ray concentration of 
the cluster. Clusters with the highest core surface brightness 
have the highest /l values. 

Like Sx, values of /l and /t do not depend strongly on 
temperature. Fig. Illl shows that the two quantities are anti- 
correlated for regular clusters, i.e. clusters with the highest 
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Figure 11. X-ray concentration parameter, /l, versus core-to- 
maximum temperature ratio, /x, for clusters at 2 = 1 and z = 0. 
Triangles are regular clusters and squares irregular clusters, with 
their median (/x, /l) co-ordinates connected by a solid line. Solid 
squares, joined by dashed lines, illustrate the median /x at fixed 
intervals of 0.1 in /l, for regular clusters. 
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X-ray concentrations have the coolest cores. As we shall see, 
these systems tend to be older (less recent merger activity), 
and thus the gas has had more time to settle down into a 
regular state. Irregular clusters tend to have low /l values 
as the subcluster boosts the overall luminosity without af- 
fecting the core luminosity. 

Examining the /l distribution alone, there is a large 
spread in values at z = 0, varying from around 0.1 — 0.7, 
but this reduces to ~ 0.1 — 0.4 at 2 = 1. As expected from 
Fig. m we see an absence of clusters with strongly-peaked 
X-ray emission at high redshift (the median /l decreases 
gradually with redshift). It is unlikely that this effect is due 
to poorer numerical resolution at higher redshift, as nearly 
all our clusters have an inner radius, Ri, that is larger than 
the physical softening radius, rmin ~ 56/i~^kpc/(l -f 2), at 
2 = and 2 = 1. Furthermore, the lack of a strong depen- 
dence of /l with temperature and the presence of cool cores 
at all redshifts suggests that numerical heating cannot be a 
major problem. 

3.6 Entropy profiles 

The combined effects of cooling and heating processes can 
be effectively probed by measuring the entropy distribu- 
tion of the ICM and comparing this with the p rediction 



of the gravitational-heatin g model (see Vpit 2005 and ref- 
erences therein). Recentlv. IVoit. Kav fc Brvan (,2005. ) com- 
pared two independent sets of gravitational-heating sim- 
ulations and found that the outer entropy profiles were 
very similar, S oc R^'^ , close to the ori ginal prediction 
from spherical accretion-shock models (e.g. lTozzi fc Normani 
120011 ). K2004a and K2004b found that cooling and feed- 
back (the same model used in this paper) only slightly 
modified the outer slope of the entropy profile; the main 
effect was an increase in the normalisation of the en- 
tropy at all radii, as suggested by o bservations (e.g. 



Ponman. Sanderson fc Finoguenovl l2003l : IPratt fc ArnaudI 
20031 ). 



The CLEF simulation allows us to study entropy pro- 
files for a much larger sample of clusters than previously. In 
Paper I, we showed that the entropy-temperature relation at 
2 = reproduced the observed behaviour of excess entropy, 
both at small (0.1ii2oo) and large (-R500) radii. In this paper, 
we focus on the scaled entropy profiles at 2 = and 2 = 1, 
shown in Fig. [I2l We scale the entropy profile of each clus- 
ter by E'^^^{z)T{^-'^^; the first factor reflects the predicted 
redshift scaling from gravitational heating, while the second 
approximately represents the scaling with temperature at 
fixed radius/overdensity, modified by non-gravitational pro- 
cesses. As will be explained below, we define Tio to be the 
projected spectroscopic-like temperature in lOkeV/fc units, 
measured between Ri and R2 (i.e. the temperature plateau, 
see subsection [33} . 

The solid curves illustrate the median and 10/90 per- 
centile profiles for regular clusters at each redshift. At both 
redshifts the profiles are close to power-law outside the core 
{R ~ (0.2 — l)ii20o); fitting a straight line to the profile 
in this region yields S cx RP'^ at both redshifts, as found 
in previous papers (K2004a,b). Note that at 7?5oo, the nor- 
malisation of the scaled profile is very similar at 2 = and 
2 = 1, thus the entropy at large radii scales with redshift as 
predicted from gravitational-heating models. 
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Figure 12. Scaled entropy versus radius, in units of -R20O1 for 
clusters at z = 1 and z = 0. The solid curves are the me- 
dian and 10/90 percentile profiles for regular clusters, the dashed 
curve is the median profile for irregular clusters and the dot- 
dashed curve for regular clusters with /l > 0.45. The verti- 
cal dashed line illustrates the force resolution of the simula- 
tion. The large shaded region encloses the mean and la stan- 
dard d eviation entropy profile for an observe d sample of clus- 
ters by IPratt. Arnaud &: Pointecouteaul l|2006h . A least-squares 
fit to the simulated data at i? > 0.2i?200 is illustrated by the 
thick solid line, fiatter than the prediction from gravitational- 
heating simulations, S oc R^-^ , shown with the lower dashed line 
llVoit. Kav fc Brvan|[2005l) . 



We also plot the median profile of irregular clusters 
(shown as the dashed curve) and at 2 = 0, clusters with the 
highest X-ray concentrations (/l > 0.45). Irregular clus- 
ters tend to have higher entropy profiles than the regu- 
lar clusters at all radii and at both redshifts. This tem- 
porary elevation in entropy reflects the shock-heating pro- 
cesses associated with the merger. Conversely, clusters with 
the highest X-ray concentrations have the lowest entropy 
proflles, reflecting the fact that they have the highest cool- 
ing rate. We also note that the profile for these systems 
resembles a broken power law, s i milar to that observed by 
iFinoguenov. Bohringer fc Zhand (|2005l ) in their REFLEX- 
DXL clusters (2 - 0.3). 

We compare the profil es at 2 = to the recent XMM- 
Newto n data studied by iPratt. Arnaud fc Pointecouteaul 
(|2006l ) (the mean plus/minus la values are shown as 
the shaded region). iPratt. Arnaud fc Pointecouteaul (|2006f ) 
use a global mean temperature, measured between 0.1 — 
0.5i?2O0, for the entropy scaling. We note, however, that 



X-ray clusters in the CLEF simulation 11 



this effectively measures their temperature plateau as they 
see no significant evidence of a de cline at large radii 
ijArnaud. Pointecouteau Sz PratJ 2005'l . The simulated (reg- 
ular clusters) and observed distributions are similar, al- 
though the simulated profile is slightly high. We note, how- 
ever, that our regular, concentrated clusters fit the observa- 
tional profile very well; the observed sample is probably bi- 
ased to systems of this type as they are intrinsically brighter 
systems and are thus easier to observe. 



3.7 Mass estimates 

Finally in this Section, we briefiy investigate the validity 
of hydrostatic equilibrium in the simulated clusters at all 
redshifts, used to estimate cluster masses from X-ray data 



Mest(< r) 



rkT{r) 



dlnp dlnT 

+ 



d In r d In r 



(7) 



where p(r) and T(r) are the 3D density and spectroscopic- 
like temperature profiles respectively. Various approxima- 
tions to equation (O have been used previously in the lit- 
erature, when little or no spatial information was available 
for the temperature distribution in clusters. Newer, high- 
quality observations with XMM-Newton and Chandra have 
overcome this problem (e.g. Arnaud. Pointecouteau fc PrattI 



l2005l : IVikhlinin et al.l l2006aV and so we assume in our study 
that the gas density and temperature profiles can be accu- 
rately recovered from the X-ray data out to -R500 • (A detailed 
study of obtaining such profiles fr om mock X-ray data is left 
to future work, but see lRasia et a l. 2006.) 

In Fig. 1131 we plot estimated to true mass ratios at iisoo 
for our clusters at z = 1 and 2 = 0. For internal consistency, 
estimated masses are those at the estimated -R500, which is 
typically 5 per cent smaller than the true -Rsoo- 

The median estimated mass is around 80 per cent of 
the true mass at both redshifts for regular clusters, with 10- 
20 per cent scatter. Irregular clusters tend to have slightly 
poorer mass estimates on average, but the scatter is also 
larger. Clusters with the highest X-ray concentrations per- 
form slightly better than the regular clusters as a whole. 

Our median mass ratio is lower than found by K2004b at 
2 = 0, who found that the average estimated mass was only 
5 per cent or so lower than the true mass at -R500, with the 
small discrepancy being due to turbulent motions (s ee also 
lEvrard. Metzler fc Navarrol Il996l : iRasia et all |2004 ). The 
reason for the difference is twofold. Firstly, K2004b used a 
mass-weighted temperature profile, where we use the more 
realistic spectroscopic-like temperature profile. This reduces 
mass esti mates by 10 percent or so, similar to what was 
found bv lRasia etall (|2006D . when assuming a low X-ray 
background in their analysis. A further 5 per cent reduc- 
tion comes from using the estimated -R500 rather than the 
true value, which also increases the scatter. While we are 
therefore not comparing true and estimated masses at the 
same radius here, we are demonstrating what the overall ef- 
fect will be on the normalisation of the mass-temperature 
relation, as will be investigated in the next section. 
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Figure 13. Ratio of estimated to true masses at -R500 versus T^i 
for clusters at 2: = 1 and z = 0. Squares are irregular clusters, 
triangles are regular clusters and filled triangles regular clusters 
with the highest X-ray concentrations, /l > 0.45. Solid (dashed) 
lines are median (10/90 percentile) ratios for regular and irregu- 
lar clusters. The dot-dashed line is the median ratio for regular 
clusters with /l > 0.45. 



4 CLUSTER SCALING RELATIONS 

We now put together the results from previous sections to 
try and understand the properties of cluster scaling relations 
in our simulation. We consider the two most important X- 
ray scaling relations in this paper: mass versus temperature 
(M — T) and luminosity versus temperature {L — T), with 
all quantities computed within J?5oo- 

Scaling relations are defined using the conventional form 



Y ^YoiX/Xoni + z f, 



(8) 



where Yo is the normalisation at A" = Xo and z = 
(for all relations, Xo = 5keV ; for the M — T relation, 
Yo = M/IO^* /i-^Mq and for the L - T relation, Yo = 
L/10''*/i~^ergs~^); a is the slope and f3 the parameter used 
to describe the redshift dependence of the normalisation. 
Scatter in the relations is measured at each redshift as 



'xiogm^ /-^$^(iog(yi/y))2 



(9) 



i.e. the r.m.s. deviation of log(y) from the mean relation, 
where Yi are individual data points. We then parameterise 
any redshift dependence of the scatter using a least-squares 

fit to Criog(y) 
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Table 1. Best-fit parameter values (and Icr errors) for scaling 
relations at z = 0. Column 1 gives the sample used in the fit; 
column 2 the best-fit normalisation; column 3 the best-fit slope 
and column 4 the logarithmic scatter. 



Sample 


'0 




o"iog{y) 


Msoo - 


Tdyn 






All 


3.44 ± 0.07 


1.38 ±0.03 


0.05 




"? Af\ -t- n ns 
o.^D m u.uo 


1 -1- n C\A 
1 .oo m u.u^ 


u.uo 




o.oo it U.ll 


i.oo it U.U/ 


U.Uo 


Low 


0.4U it U.ii 


i.oo ± U.Uo 


U.Uo 


Msoo - 


Tgas 






All 


4.08 it 0.06 


1.69 ±0.03 


0.03 


Reg 


4.05 ±0.06 


1.67 ±0.03 


0.03 


High 


3.95 ±0.12 


1.59 ±0.06 


0.03 


Low 


4.07 ±0.08 


1.69 ±0.03 


0.03 


Msoo — 








All 


5.51 ±0.32 


1.81 ±0.08 


0.08 


Reg 


4.93 ±0.29 


1.69 ±0.08 


0.08 


High 


7.03 ± 0.78 


1.96 ±0.16 


0.06 


Low 


4.37 ±0.21 


1.61 ±0.07 


0.06 


Msoo - 


T^SO 

-'si 






All 


4.47 ±0.19 


1.76 ±0.07 


0.08 


Reg 


4.02 ±0.12 


1.64 ±0.05 


0.05 


High 


3.94 ±0.13 


1.62 ±0.07 


0.03 


Low 


4.08 ±0.17 


1.66 ±0.06 


0.05 



iw-est 
™500 


T^SO 
-"si 






All 


3.28 ±0.15 


1.64 ±0.08 


0.09 


Reg 


3.11 ±0.11 


1.57 ±0.06 


0.06 


High 


3.22 ±0.15 


1.56 ±0.10 


0.04 


Low 


2.98 ±0.15 


1.53 ±0.08 


0.06 



L\,oi - 

All 

Reg 

High 

Low 



6.0 ± 1.1 
5.8 ±1.3 
19.0 ±6.5 
3.6 ±0.5 



3.08 ±0.26 
3.04 ±0.30 
3.61 ±0.45 
2.81 ±0.19 



rSO 

AU 
Reg 
High 
Low 



7^50 
■ -'si 



4.3 ±0.3 
3.9 ±0.3 
4.7 ±0.3 
3.1 ±0.3 



3.45 ±0.13 
3.37 ±0.12 
3.16 ±0.12 
3.13 ±0.14 



(o-log(y)) {z) = (To + (71 log(l + z). 



0.27 
0.27 
0.16 
0.16 



0.14 
0.12 
0.05 
0.11 



(10) 



4.1 Results at 2; = 

Our results for z = G clusters are summarised in Table [1] 
Column 1 lists the sample used when fitting the data. Here, 
we consider all 95 clusters in our temperature-limited sam- 
ple (labelled Alt), the 83 regular clusters (i.e. those with 
S'x ^ 0.1; labelled Reg), the 23 regular clusters with the 
most prominent core emission (/l > 0.45); denoted High) 
and the 60 remaining regular clusters (denoted Low). Col- 
umn 2 lists the best-fit normalisation, lo; column 3 the slope 
of the relation, a; and column 4 the scatter in the relation, 
f"iog(y)- We now discuss each relation in turn. 
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Figure 14. Estimated scaled-mass at R$oo versus spectroscopic- 
like temperature outside the core at z = 0. Squares are irregular 
clusters, triangles are regular clusters and filled triangles regular 
clusters with the highest X-ray concentrations, /l > 0.45. The 
solid line in each pa nel is a best fit to regular clusters. Cr osses 
are data-points from | Arnaud. Pointecouteau &: Pratj |200^ and 
IVikhlinin et al.l l|2006a^ . 
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We first study the M ~ T relation at z — 0. In Pa- 
per I we presented results for the hot gas mass-weighted 
temperature within i?2500, Tgas, and showed that the re- 
lation was in good agreement w ith the Chandra results of 
lAUen. Schmidt fc FabianI (|200ll ). Here we discuss the M- 
T relation at 7?5oo as we expect it to be less suscepti- 
ble to cooling and heating effects associated with the clus- 
ter core. While measuring the relation at J?soo is obser- 
vationally challenging, even with XMM-Newton and Chan- 
dra, recent attempts have been performed for a small sam- 
ple of clusters with a reasonable ra nge in temperature 
llArnaud, Pointecouteau fc Prat3 |2005| . hereafter APP05; 
IVikhhnin et al.l bOOGaT hereafter VKF06). We will eventu- 
ally compare our results at z = to these observations, 
but first study how our definition of temperature and mass 
affects the details of the relation. 

We initially consider the relation between the true to- 
tal mass of a cluster, Msoo, and its dynamical temperature, 
Tdyn, where the latter was defined in equation ((3}. This re- 
lation should most faithfully represent the scaling expected 
from gravitational-heating models (a = 1.5) but as listed 
in Table [1] the measured slop e is slightly shallower than 
this (a ~ 1.4). As discussed in iMuanwong. Kav fc Thomas! 
(2006), the deviation in slope is consistent with the varia- 
tion in halo concentration with cluster mass (i.e. even the 
dark matter is not perfectly self-similar). Note that the sub- 
samples give almost identical results to the overall sample, 
although the High sub-sample exhibits less scatter. 

We next consider the hot gas mass-weighted tempera- 
ture, Tgas- The slope of the relation steepens to a ~ 1.7; as 



discussed in Section [3.11 this is due to the combined effects 
of heating and cooling. Strikingly, the scatter in this relation 
is very small ((Tiog{T) ~ 0.03). Again, no significant change 
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Figure 15. Bolometric luminosity versus spectroscopic-like temperature, for all emission (left panel) and emission outside the core 
(right panel). Squares are irregular clusters, triangles are regular clusters and filled triangles regular clusters with the highest X-ray 
conc entrations, > 0. 45. The solid line in each pan el is a best fit to regular clusters. Data points with error bars are observed values 
from lMarkevitchI lll998h and lArnaud fc EvrardI ^9^. 



in the relation is observed when the cluster sub-samples are 
considered. 

When the X-ray temperature, Tsi, is used, both the nor- 
malisation and scatter increase, with the irregular and High 
clusters lying above the mean relation (i.e. they are colder 
than average). This is because cool, dense gas in the core 
and in substructures throughout the cluster (see Fig. [2]) is 
weighted more heavily than before, and there is a large vari- 
atio n in the cool gas distribution from cluster to cluster (se e 
also iMuanwong. Kay fc ThomasllioO^ : lO'Hara et al.ll2006h . 
As can be seen in Fig. 1111 the irregular and the High clusters 
have the lowest /t values. 

We also present results for the spectroscopic-like tem- 
perature when particles from within the inner 50/i"^kpc 
core are excluded (denoted Tj"), which reduces the scatter 
in the Reg clusters from 0.08 to 0.05. The High and Low 
relations are now consistent with the overall Reg relation, 
although the irregular clusters still lie above the relation as 
a second cool core is still present. 

Finally, we replace the actual mass with the mass es- 
timated under the assumption of hydrostatic equilibrium 
(denoted M|oo), as defined in Section [331 Fig. [U illus- 
trates the result, in comparison to the data from APP05 
and VKF06. The relation for the Reg subsample provides 
the closest match to the observational data. The main effect 
of using the estimated mass is to reduce the normalisation by 
~ 20 per cent. Although the two observational samples are 
similar, our Reg relation is closest to the best-fit results of 
VKF06; the slope and scatter are almost identical (VKF06 
find a = 1.58 and (Tiog(A/) ~ 0.06) and the normalisation 
differs by 10 per cent or so (VKF06 find Yo = 2.89 ± 0.15). 
Given the variations between parameters considered in this 
study, this is quite a good match, but serves to point out 
that a precision measurement of the M — T relation is non- 
trivial and must include several physi cal effects; our pre sent 
study is by no means exhaustive (see lRasia et al.ll2006l ). 



4.1.2 L - T relation 

To study the luminosity-temperature relation, we compute 
bolometric luminosities, Lboi, for all emission within iisoo 
(where more than 90 per cent of the cluster emission comes 
from). We also compute luminosities outside the core (de- 
noted I/boi), again by excluding all hot gas particles from 
within 50 h~^kpc from the cluster's centre. 

Fig. illustrates luminosity-temperature relations at 
2 = 0. In the left panel we show results for total luminosities 
and spectroscopic-like temperatures, and in the right panel, 
for luminosities and temperatures outside the 50^~^kpc 
core. Best-fit parameters for the various cluster samples at 
z — are also given in Table [1] 

When all emission is included, the L — T relation at 
z — has a large amount of scatter. Comparing the relation 
for regular clusters with high /l values to those with low 
/l values, we see that the two subpopulations are widely 
separated in the L — T plane. The scatter thus reflects the 
strength of the cor e emission, as shown observationally by 
iFabian et al.l (|l994|). We discuss this further in Section [S] 

When the 50/i"^kpc core emission is excised, the scatter 
in the relation reduces substantially, from 0.27 to 0.14, with 
all samples then having very similar properties. We also note 
that irregular clusters do n ot lie systematically off the L — T 
relation, in agreement with lRowlev. Thomas fc Kavl (|2004 ). 
who analysed a simulation with radiative cooling but no 
feedback. 

We compar e our excised-core resul ts with the obser- 
vation al data of iMarkevitchI (| 19981 ) and lArnaud fc EvrardI 
(| 19991 ): the former also excised emission from the inner 
50/i~^kpc and the latter selected non-cooling-flow clusters. 
Although our clusters do not cover the same dynamic range 
as the observations, we note that our L—T relation has a nor- 
malisation that is too high (see also Paper I). This suggests 
that that cluster temperatures in general are too low (note 
that higher temperatures may not signiflcantly affect the 
normalisation of the M — T relation, as the estimated mass 
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depends linearly on T) . The slope of the relation for all clus- 
ters, a — 3.5, is steeper than the observations (a ^ 2.6—2.9). 
As stated in Paper I, the slope varies systematically with 
temperature, such that higher-mass clusters have lower val- 
ues. The lack of hot clusters in our sample biases our result 
to higher values. Significantly larger volumes are still re- 
quired to capture the rich clusters, to get a more accurate 
(average) slope for the cluster population. 



4.2 Evolution of scaling relations 

We now study how the M — T and L — T relations evolve 
with redshift. We first measure the slope, normalisation and 
scatter of the relations at each redshift between z — and 
z = 1. The gravitational-heating model predicts the slope 
to be constant with redshift. For the M — T relations this 
is generally true; although the variation can be quite noisy, 
there is no evidence for a systematic change in the slope, a, 
with redshift (e.g. see the top-left panel in Fig.[Tn]for how the 
slope changes with redshift in the Mfoo — T^f relation) . For 
the L — T relation (all emission), the slope increases with 
redshift when all clusters are considered. This is because 
the few hottest clusters have anomalously high temperatures 
for their luminosity at low redshift, causing a decrease in 
slope since they carry a lot of weight. At higher redshift the 
effect diminishes as the clusters move back towards the mean 
relation. We circumvent this problem by restricting our fit 
to the L — T relation to clusters with 2 < kT < 5 keV at each 
redshift; as can be seen in the top-right panel of Fig. 1161 the 
slope of the I/boi— Tsi relation is now approximately constant. 
For all relations, we fix a to its median value between z = 
0-1. 

With Q determined, we then fit equation ^ to the nor- 
malisation data to determine Yb and /3. (Note this may cause 
Yo to change slightly from the exact z = Q value.) The scat- 
ter is also determined at each redshift (equation [5| and fit 
with equation (|10p . 

Table [2] gives best-fit parameters for our generalised 
scaling relations when applied to all clusters at each red- 
shift. For the E[z)M — T relations, we see a lack of evolution 
relative to the simple scalings predicted from gravitational 
heating, with |/3| < 0.15. The scatter also changes very little 
with redshift, with \a\ \ < 0.1 in all cases. This lack of evolu- 
tion in normalisation and scatter is illustrated more clearly 
for the E[z)Mifo - Tf relation in the left panels of Fig.fTUl 

The evolution of the E~^{z)L — T relation is also pre- 
sented in Fig. [16] (see also Table [l]). Contrary to the M - T 
relation, this relation evolves negatively with redshift, with 
/3 ~ — 1. Note the amount of evolution at jz = 1 is compa- 
rable to the intrinsic scatter in the relation at 2: = 0. What 
is striking from the figure, however, is the evolution of the 
scatter with redshift: o"iog(L) at z = 1 is almost a factor of 3 
lower than at z = 0. As was found in subsection [331 the dis- 
persion in X-ray concentration decreases with redshift, such 
that at high redshift, clusters with strong cooling cores are 
absent. This is refiected here as a reduction in the scatter 
of the L — T relation. When the core is excised, the scatter 
is reduced at all redshifts and also evolves less. The normal- 
isation also evolves less with redshift, demonstrating that 
some (but not all) of the deviation from the gravitational- 
heating case is due to processes occurring within the inner 
core. Furthermore, since we know that the E{z)M5oo — T^i 



relation evolves very weakly with redshift, negative evolu- 
tion m the E'^{z) 

Lhoi ~ '^ai' relation is almost entirely due 
to a deficit in luminosity, again as seen in the entropy and 

surface-brightness profiles. 

A similar study was perforr ned by Ettori et al. J 2004 ) , 
using the same simulation as iBorgaiii et al.l 1 2004 ). Al- 
though they used a difi'erent model for star formation 
and feedback than used here, they obtained very sim- 
ilar results for the evolution of the E(z)M — T and 
E^^{z)L — T relations; using our notation, they found 
P — —0.2 and /3 = —0.8 re s pectiv ely. On the other hand, 
iMuanwong. Kav fc Thomas! (|2006l ) compared a simulation 
similar to (but smaller than) the CLEF simulation, with 
a simulation with radiative cooling only and with a sim- 
ulation with cooling and preheating. They found that the 
evolution of the L — T relation varied enormously between 
the models. Their conclusion was that the amount of evolu- 
tion depended on the nature of non-gravitational processes. 
We can thus conclude, at this point, that no general con- 
sensus has emerged from numerical simulations as to what 
the expected evolution of cluster scaling relations will be, 
once sufficiently-large samples of high-redshift clusters ex- 
ist. Of vital importance, from the simulation side, will be 
to produce cluster catalogues that are well matched to the 
observations; in particular, the deficit of high-temperature 
systems in most studies to date needs to be addressed. 



5 DISCUSSION 

Perhaps the most interesting result in this paper is that 
our simulation predicts a large scatter in the luminosity- 
temperature relation at low redshift, as observed, however 
this scatter decreases with redshift due to the lack of sys- 
tems with high X-ray concentrations at 2: ~ 1. Here, we 
discuss this issue in more detail and investigate further the 
differences between clusters with high and low X-ray con- 
centrations, and cool and warm cores. 



5.1 Mass deposition rates 

Observed samples of (generally low-redshift) clusters are his- 
toric ally split into cooli ng-flow and non- cooling- flow systems 
(e.g. iFabian et ai]|l994h , with the former having higher mass 
deposition rates, usually estimated from their core luminos- 
ity and temperature 



2 /imni 
5 kT 



(11) 



X-ray spectroscopy of cluster cores has revealed that sig- 
nificantly less gas in high Mx clusters is actually cooling 
down to temperatures significantly below the mean temper- 
ature of the cluster. This lack of cold gas is likely attributed 
to intermittent hea ting from a c entral active galactic nu- 
cleus (AGN; e.g. see lFabianir2003l for a recent review) . How- 
ever, it is important to understand the origin of the large 
spread in Mx within the cluster population, as it also ex- 
plains m uch of the scatter in the luminosity-temperature 
relation (|Fabian et al.lll994 ). 

We have measured Mx for our clusters (within a fixed 
physical radius of rcorc ~ 50 /i~^kpc) and, as expected, 
found that it is strongly correlated with /l, ranging from 
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Figure 16. Redshift dependence of the slope, normalisation and scatter of the M|qq — T^^ (left panel) and Lbol — T^i (right panel) 
scaling relations. The band in the top panel illustrates the mean slope plus the standard deviation at each redshift, all normalised to 
the median value over all redshifts between 2 = — 1. The band in the middle panel illustrates the mean and standard deviation of the 
normalisation at each redshift (assuming a fixed slope, at the median value). The best-fit straight line to the mean data is also plotted as 
a dashed line and is used to normalise the data at 2 = 0. The band in the bottom panel illustrates the logarithmic scatter in the scaling 
relations at each redshift (values are reflected about the x-axis to give an idea of the full size of scatter in the relation), with the best-fit 
straight line given by the dashed line. 



Table 2. Best-fit parameter values (and la errors) for evolution of scaling relations from 2 = — 1. Column 1 gives the median slope 
used for the fit; columns 2 &; 3 the best-fit normalisation and evolution parameters; and columns 4 & 5 the best-fit scatter parameters. 



(a) Yo 


13 


CTO 


CTl 


E{z)M500 - Tdyn 

1.40 3.48 ±0.02 


-0.07 ±0.01 


0.05 ± 0.002 


-0.02 ±0.01 


E{z)M500 - Tgas 

1.68 4.08 ±0.02 


-0.12 ± 0.01 


0.04 ± 0.004 


0.03 ± 0.02 


E{z)M5oo - %i 
1.77 5.21 ±0.04 


-0.13 ± 0.02 


0.08 ± 0.003 


-0.05 ±0.01 


E{z)M5oo - 

1.67 4.24 ±0.03 


0.04 ±0.02 


0.07 ±0.003 


-0.03 ±0.02 


E{z)MI-l - 

1.56 3.17 ±0.04 


-0.07 ±0.03 


0.08 ±0.02 


0.08 ± 0.10 


E-i(^)iboi -TbI 
3.36 7.36 ± 0.09 


-0.98 ± 0.03 


0.27 ±0.003 


-0.59 ±0.02 


^ bol si 

3.41 4.53 ± 0.07 


-0.61 ±0.04 


0.14 ±0.004 


-0.23 ± 0.02 
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Figure 17. Inferred mass deposition rate from X-ray emission, 
Mx, versus X-ray concentration for regular (triangles) and irreg- 
ular (squares) clusters at 2 = 1 and 2 = 0. Filled triangles are 
regular clusters with /l > 0.45. 



~ 1 - 900AfQyr"i (Fig. [TT]). Clusters with the highest 
concentrations, /l > 0.45, are regular and typically have 
Afx > lOOA'/eyr^^. These clusters could be called strong 
cooling-flow systems as they most resemble the observational 
samples of the same name; again note the absence of these 
objects at z = 1. 

While the median /l decreases with redshift, the me- 
dian Afx stays approximately constant. The lack of strong 
cooling-flow clusters at high redshift is offset by the increase 
in Mx for individual systems, due to the ratio, rcoro/iisoo, 
being typically larger at higher redshift, thus capturing more 
of the cluster's luminosity. Averaging over all redshifts, we 
found that (Mx) = 35 ± 5M0yr"\ 

Here, we do not attempt to address the issue of how 
much gas is actually cooling down within our cluster cores. 
As discussed in K2004b, our simulations currently lack the 
number of particles to accurately follow the inward flow of 
the gas all the way down to low temperature. However, as 
we will demonstrate below, we find that the large range in 
X-ray concentration/cooling-flow strength exhibited by our 
clusters at low redshift is strongly dependent on the clus- 
ter's larger-scale environment, i.e. whether it experienced a 
late-time major merger or not. So while the dynamics of a 
cooling core within a given cluster may not be accurate, and 
requires further investigation, our main (statistical) conclu- 
sions should hold as the simulation has accurately followed 
the merger histories of the cluster population. 
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Figure 18. Redshift when each cluster last had a warm core 
versus the nearest redshift when it was irregular. Clusters that 
were regular (irregular) at 2warm are shown as triangles (squares). 
The solid line is 2warm — ^irr • 



5.2 Cooling flows, cool cores and dynamical state 

Besides their high core luminosity, cooling-flow clusters have 
traditionally assumed to be dynamically-relaxed systems 
hosting a cool core. Conversely, non-cooling-flow clusters 
with low core luminosities are thought to host isother- 
mal/warm cores and be dynami cally disturbed. This vi ew- 
point was recently challenged bv lMcCarthv et al.l l|2004h as 
being overly-simplistic, as observations of both cooling-flow 
clusters with disturbed morphologies (e.g. Perseus) and non- 
cooling-flow clusters (e.g. 3C 129) with relaxed morphologies 
exist. Our simulation lends some support to their argument, 
as Fig. 1171 shows. At z = 0, irregular clusters are found to 
have a large range in /l (or A^x), with one irregular cluster 
{Sx = 0.14) having Afx = 358AfQyr"^. Conversely, reg- 
ular clusters can also have very low X-ray concentrations 
(Afx < lOAfeyr"^). However, statistically, the average reg- 
ular cluster has a higher X-ray concentration than an ir- 
regular cluster. This is because the X-ray concentration is 
related to the dynamical history of the cluster, as we will 
show below. 

We showed in subsection 13.51 that /l is anti-correlated 
with the strength of the cool core, /t, as measured from the 
projected temperature profile; clusters with the coolest cores 
have more concentrated X-ray emission. However, nearly 
all of our clusters, regular and irregular, have cool cores 
(/t < 1; Figs. llOlfc lll|l . This is in agreement with previous 
simulation work wh ere the gas was allowe d to cool radia- 
tivcly (c. g.lMotl et al. 2004; Rowley, Tho mas fc Kwll2004l : 



IPooleet al.[ 20061 ) . where it was found that cool cores are 



very hard to disrupt by mergers. 

Warm (or non-cool; /t 1) cores exist but are rare 
in our simulation. Given the number of outputs available, 
only one quarter of the clusters were found to host a warm 
core since z — 1, lasting at most around 1 Gyr. Interest- 
ingly, clusters with warm cores nearly always appear reg- 
ular, even though the generation of a warm core appears 
linked to the merger process. This is shown in Fig. 1181 
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Figure 19. Distribution of /x values at 2 = (solid 
histogram), compared to the obs ervational sample of 
[Sanderson. Ponman fc O'SuUivanl 1I2OO6I '). 



where we see a clear correlation between the redshift when 
a cluster last had a warm core, Zwarm, against the near- 
est redshift when it was irregular (Sx > 0.1), zin. It 
is unclear whether the cores are heated solely from the 
gravitational interaction of the merger, or a contribution 
comes from the feedback, which could also be triggered 
by a merger. Nevertheless, the paucity of warm cores is 
at odds with the observational data at low red shift. For 
example, ISanderson. Ponman fc O'SullivanI (|2006l ) recently 
studied a flux-limited sample of 20 clusters observed with 
Chandra, and found only half of them to contain cool 
cores, even though the core gas in the warm core clus- 
ters have cooling times significantly shorter than a Hub- 
ble time. The discrepancy is illustrated clearly in Fig. 1191 
where we compare the distribution of /x values found in 
t he CLEF simulation at z = with the o bservational data 
of lSanderson. Ponman fc O'SullivanI (|2006l ) . Although based 
on a limited sample, the observations suggest there exists a 
bimodal distribution, not present in the simulation. This 
suggests that our simulation is still missing a heating mech- 
anism that could produce a larger fraction of warm cores, 
which again could be linked to AGN activity. 



5.3 Scatter in the L— T relation 

We now examine why there is a large scatter in the L — T 
relation at low redshift. Classically, it is thought that the 
scatter is related to the dynamical histories of clusters. In 
particular, clusters with the strongest cooling flows (which 
lie above the mean L — T relation) are believed to be in that 
state because they have not endured a major merger in the 
recent past. Our simulation supports this picture, as will be 
demonstrated in the following two figures. Firstly, Fig. l2Ul ex- 
plicitly shows that the scatter in the L — T relation is tightly 
correlated with the X-ray concentration (or mass deposition 
rate) of a cluster. For regular systems, clusters with higher 
X-ray concentrations lie above the mean relation, and those 
with low X-ray concentrations below. Irregular clusters lie 



Figure 20. Offset in luminosity from the mean L — T relation ver- 
sus X-ray concentration for clusters at 2 = 0. Triangles are regular 
clusters and squares irregular clusters. Solid triangles are clusters 
with the highest X-ray concentrations/X-ray-inferred mass depo- 
sition rates. 
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Figure 21. X-ray concentration at z = versus the lowest 
redshift when the cluster was irregular. Only clusters in our 
temperature-selected sample at all available redshifts, z < 2;^ 
are plotted. Triangles are regular clusters today and squares ir- 
regular clusters {z^^^ = 0). Asterisks are clusters with z^^^ > 1, i.e. 
they did not experience a major merger between now and 2 = 1. 
The dashed line is a best-fit relation to the regular clusters with 

2irr < 1. 



off this correlation because /l decreases due to the presence 
of a second object (which also boosts the luminosity). 

Secondly, in Fig. [21] we plot X-ray concentration at 
2 = versus the lowest redshift when each cluster expe- 
rienced a major merger. Clusters which are not present in 
our temperature-selected samples at all redshifts, z < Ziir 
are not plotted. Clearly there is a strongly-positive corre- 
lation, demonstrating that the most concentrated systems 
did not experience a major merger in the recent past (the 
asterisks are those clusters with Zin > 1). 
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An alternati ve mechanism fo r gene rating the scatter 
was proposed by [McCarthy et ahl l|2004 ). who used semi- 
analytic models of clusters with preheating and cooling (but 
the effects of accretion and merging of haloes were not in- 
cluded). They suggested that the position of a cluster on the 
luminosity-temperature relation was related to the level of 
preheating it received: clusters that experienced higher lev- 
els of preheating correspond to non-cooling-flow clusters (i.e. 
low X-ray concentratio ns, here) and vice-versa. Similarly to 
iMcCarthv erall l|2004l ). we tested whether the amount of 
feedback is correlated to the strength of the cooling core. 
Such an effect should be seen through a trend of stellar 
mass fraction with /l, as our feedback model injects energy 
approximately in proportion to the star-formation rate. No 
trend is seen in our simulation, i.e. stellar mass fractions are 
similar between clusters with low and high X-ray concentra- 
tions. 

It is clear, therefore, that the strong cooling-flow popu- 
lation exists in our model at low redshift because of a lack of 
major merger activity in such systems &t z < 1. The absence 
of strong cooling-flow systems at higher redshift, responsi- 
ble for the decrease in the L — T scatter, can therefore be 
attributed to the increase in the merger rate with redshift. 

The absence of strong coohng-flow clusters at high red- 
shift in our model has important implications for cluster 
cosmology. Large samples of X-ray clusters at high redshift 
are still in their infancy, although will start to become avail- 
able over the ne xt few years, such a s from the XMM-Newton 
Cluster Survey (|Romer et al.ll200ll '). 

If our prediction is correct, it will have both positive and 
negative implications for cosmology. On the positive side, 
the smaller scatter will allow for a simpler survey selection 
function, with incompleteness effects being less of a problem. 
On the negative side, there is a lack of very luminous objects, 
so the number of high-redshift clusters above a given flux 
limit will be considerably less, reducing the overall power 
for specific surveys to constrain cosmological parameters. 
Interestingly, first observational results see m to support the 
lack of cooling-flow systems at high redshift (jVikhlinin et al.l 
l2006bl ). 

Another interesting point that our result throws up, is 
whether strong cooling-flow clusters would exist in a uni- 
verse with Qni = 1? In such a model, the merger rate would 
be expected to change very little with redshift, so clusters 
today may not have had the time to establish a strong cool 
core. In other words, the strongest cooling-flow clusters only 
exist because of the freeze-out of structure formation in a 
universe with sub-critical matter density. 



6 CONCLUSIONS 

In this paper, we presented the cluster population that forms 
within the CLEF simulation, an A''-body/hydrodynamics 
simulation of the ACDM cosmology, with radiative cooling 
and energy feedback from galaxies. Our cluster sample, with 
nearly one hundred kT > 2 keV objects at z = and sixty 
at 2 = 1, is one of the largest drawn from a single simu- 
lation. In this paper, we studied the demographics of the 
cluster population out to z = 1, focusing on the effects of 
dynamical activity and the strength of cooling cores, and 
how the X-ray properties of clusters depend on them. The 



Sunyaev-Zel'dovich properties of the clusters may be found 
in a companion paper (da Silva et al., in preparation). Our 
main conclusions are as follows; 

• We quantified the amount of dynamical activity (ma- 
jor mergers) within the cluster population, using a simple 
projected substructure statistic, based on the observable X- 
ray surface-brightness distribution. While there is no signif- 
icant dependence of this quantity, Sx, with cluster temper- 
ature, it does increase with redshift. The fraction of irregu- 
lar, Sx > O.I clusters, shown to be merging systems in the 
surface-brightness maps, increases from around 10 per cent 
at z = to 20 per cent at z = 1 , thus constituting a minority 
population at all redshifts. 

• The projected ICM temperature profile of regular clus- 
ters has a generic shape at low and high redshift, decreasing 
in the centre (due to radiative cooling) and beyond 0.2i?5oo, 
due to the intrinsic shape of the gravitational potential. Ir- 
regular clusters have flatter profiles at large radii due to the 
presence of a second object which compresses and heats the 
gas. The shape of the regular cluster profile at z = is in 
go od agreement w i th the recent study of cool core clusters 
bv lVikhhnin et all l|2005l ). 

• To quantify the core properties of our clusters, we de- 
fined two simple (and observationally-measurable) struc- 
ture parameters, /t, which measures the core to maxi- 
mum temperature ratio, and /l which measures the frac- 
tion of emission from within the core (the X-ray concen- 
tration of the cluster). We found that the vast majority 
of clusters contain cool cores (/t < 1) at all redshifts. 
This is at odds with the observational data, at least at 
low redshift, where only half of clusters contain cool cores 
l|Sanderson. Ponman fc O'SuUivanl l2006l ). The X-ray con- 
centration, /l, is anti-correlated with /t. The dispersion in 
/l is large at z = 0, but decreases with redshift due to the 
absence of clusters with the highest values (i.e. the strongest 
cooling cores). 

• The scaled entropy profile has an outer logarithmic 
slope of 0.9 and decreases all the way into the centre, with 
no evidence of a flattened core. The ratio of the normal- 
isation at large radii, for clusters at 2 = 1 and z = 0, 
is similar to that expected from the gravitational-heating 
model (^(r) cx E~*/^{z)), but the z = 1 clusters have 
higher central entropy than at 2 = 0. Irregular clusters 
have higher entropy profiles and regular clusters with strong 
cooling cores have lower entropy profiles. The profile at 
2 = (in particular for the strong cooling core clusters) 
is in good agreement with the recent observational data of 
IPratt. Arnaud fc Pointecouteaul l|2006l '). 

• Mass estimates of X-ray clusters, based on the hydro- 
static equilibrium equation, are around 20 per cent lower 
than the true masses, even when spatial density and tem- 
perature in f ormat ion of the ICM is known. As found by 
iRasia et al.l (|2006| ). the reasons for the discrepancy are X-ray 
temperature bias to low entropy gas and incomplete ther- 
malisation of the gas. 

• The estimated mass versus spectroscopic-like temper- 
ature relation at 2; = is only ~ 10 per cent higher than 
the observed relation for 7?5oo. Splitting the regular cluster 
sample into those with weak and strong cooling cores makes 
little difference to the properties of the relation, when the 
temperature is measured outside the core. Thus, details of 
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the mass-temperature relation should be insensitive to the 
cluster selection procedure. 

• The mass-temperature relation evolves similarly to the 
gravitational-heating model prediction, M{T) oc E~^{z). 
The scatter, A(logM) ~ 0.08, evolves very little with red- 
shift. 

• The luminosity-temperature relation has a large degree 
of scatter at z = 0, reflecting the large dispersion in X-ray 
concentration of the clusters. Excising the core emission re- 
duces the scatter considerably, although leads this to a re- 
lation that still has a higher normalisation than observed. 
Irregular clusters are not systematically offset from the main 
relation. The luminosity-temperature relation evolves neg- 
atively with rcdshift, contrary to the gravitational-heating 
expectation, where L{T) oc E{z). Excising the core reduces 
this negative evolution, with almost self-similar evolution at 
very low redshift. 

• The scatter in the luminosity-temperature relation de- 
creases strongly with redshift, again due to the lack of strong 
cooling core clusters at high redshift. There is a positive cor- 
relation between the X-ray concentration of the cluster and 
the redshift when it last had a major merger, but apparently 
not between the X-ray concentration and the level of feed- 
back experienced by the cluster. Thus, our results indicate 
that the formation of a cooling-flow population of clusters at 
low redshift is tied to the slow down in dynamical activity in 
the ACDM model, allowing clusters in quieter environments 
to develop a strong cooling core. 

Our simulation is one of the first of a new generation 
that is able to follow a substantial rmmbcr of objects with 
reasonable resolution, while attempting to include the vi- 
tal physical processes that alter the gravitationally-heated 
structure of the ICM: radiative cooling, star formation and 
feedback. While our particular model can reproduce many 
observed characteristic features of the cluster population, 
particularly those with cool cores, we acknowledge that it 
has its shortcomings. For example, it fails to completely 
quench the overcooling of baryons into stars, it does not 
predict enough clusters with warm cores, and it docs not 
match the L — T normalisation in detail (being too high). 

All these problems point to the need for an even more 
efficient heating mechanism that reduces further the amount 
of cool gas in the clusters, without destroying the already 
good agreement in cool core clusters. It may be possible 
that the problems could be overcome by fine tuning the two 
feedback model parameters. However, it is desirable to incor- 
porate a more realistic physical model for feedback, that is 
able to treat separately the effects from stars and black holes 
(in our current model, the heating rate directly follows the 
star- formation rate). The wealth of high-quality X-ray data 
that is becoming available will undoubtedly help constrain 
the feedback physics further, and thus allow more realistic 
cluster models to be constructed. 
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